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Abstract 

We predict new generic types of vorticity-carrying soliton complexes in 
a class of physical systems including an attractive Bose- Einstein conden- 
sate in a square optical lattice (OL) and photonic lattices in photorefrac- 
tive media. The patterns include ring-shaped higher-order vortex solitons 
and supervortices. Stability diagrams for these patterns, based on di- 
rect simulations, are presented. The vortex ring solitons are stable if the 
phase difference A(j> between adjacent solitons in the ring is larger than 
7r/2, while the supervortices are stable in the opposite case, A(j> < n/2. A 
qualitative explanation to the stability is given. 

Recent experimental observations of two-dimensional (2D) solitons and 
localized vortices [5] in photonic lattices in self-focusing photorefractive crystals 
(PRCs), and prediction of similar structures in the 2D Gross-Pitaevskii equation 
(GPE) for a self-attractive Bose-Einstein condensate (BEC) loaded into the 
square optical lattice (OL) 0] , have drawn interest to solitons of this type, 
including vorticity-carrying ones. Vortex solitons of the gap type have also been 
predicted in the repulsive GPE with the square OL |5lll0|. including higher-order 
ones, with the vorticity S = 2 |1(J) . However, the results which were reported 
thus far for the self-focusing media, both saturable and cubic, pertain solely to 
the vortices with S = 1 (including strongly deformed ones |ll)b In this work, 
we present stable vortex solitons with higher values of S and more complex 
vortex patterns, which suggests a possibility to create them in the experiment. 
We report the results for both BEC and PRC models, which attests to their 
generic character. It is relevant to mention that, in discrete models, which 
may be considered as a limiting case corresponding to a very strong OL, with 



attraction and repulsion alike, both 2D and 3D vortex solitons with S = 2 are 
unstable, but the ones with 5 = 3 have their stability region J3| • 

In a moderately strong OL, we find stable vortex solitons with higher vortic- 
ities S and their counterparts in the form of quadrupole solitons (in the discrete 
model, the latter ones may be stable too [T2]l. The vortices with large S fea- 
ture a ring shape, resembling vorticity-carrying "soliton necklaces" , that were 
recently studied in diverse spatially uniform J3j models, but, on the contrary 
to those models, the ring-shaped vortices on the square lattice are true sta- 
tionary states, rather than slowly disintegrating quasi-patterns in the free space 
(recently, zero-vorticity necklaces and other soliton complexes were predicted 
and observed in Kerr and PRC media with a photonic lattice We also 

construct another generic type of stable structures that we call "supervortices" 
|15|. i.e., ring-shaped (or densely packed) arrays of individual compact vortices 
with 8 = 1, arranged so that the ring itself carries global vorticity (s pertains 
to the intrinsic vorticity of an individual vortex in the array, to distinguish it 
from S referring to the global vorticity imprinted into the array) . Although the 
ring-shaped supervortices may formally resemble dark-soliton vortex necklaces 
in external traps |16j , they have never been considered before, as the individual 
solitary-wave vortex in unstable without the OL. 

In the normalized form, the two-dimensional (2D) GPE with the negative 
atomic scattering length and square-OL potential is [3] 

iipt + (l/2)V 2 i/> + \^\ 2 ^P + £ [cos(2x) + cos(2y)] tp = 0, (1) 

where ip(x,y,t) is the BEC wave function, the OL period is scaled to be ir, 
and e is the OL strength, in units of the respective recoil energy. The vertical 
coordinate is excluded, assuming tight confinement in the third direction, which 
is readily provided by light sheets 0. As is known (see, e.g., Ref. [H]), the 
underlying three-dimensional GPE with a negative scattering length of atomic 
collisions, a < 0, can be reduced to its 2D counterpart JJJ, provided that the 
transverse quantum pressure is much stronger than the self-attraction (which 
precludes the collapse). In other words, the energy of the atom with the mass 
m in the ground state of the vertical trap with a size a z , E { Z) ~ h 2 /(2ma 2 z ), 
must be much larger than the contribution of the attraction between atoms 
to the atomic chemical potential, A/i ~ —lnh 1 (o/m)n, where n is the three- 
dimensional density of atoms. Thus, the 2D GPE equation for the self-attractive 
condensate is relevant under the condition n <C (47ra^a) , which can be readily 
met in the experiment (in particular, a can be very small in 7 Li |7j). In fact, 
this condition can be relaxed due to the stabilizing effect of the OL potential in 
Eq. JTJ. 

Equation Q does not include the external trapping (parabolic) potential, as 
we aim to look for patterns whose stability does not depend on the presence of 
the trap. An equation governing the spatial evolution of the amplitude ip(x, y, z) 
of the probe beam in the PRC with a square photonic lattice of strength I differs 
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from Eq. by the presence of saturation ^IS], 

**> + (1/2)V ^ + l + / [cos(^) + cos(^)] 2 + |^p (2) 

(strictly speaking, PRCs feature anisotropy, but in crystals used in work 1 
the anisotropy is very weak and may be neglected). Stationary solutions with 
chemical potential ^ are looked for as i/j(x, y, t) = e~' lfJ ' ti S(x, y), where the func- 
tion \E' is, generally, complex (in the PRC model, — fx is the propagation con- 
stant). In this notation, the solutions depend on two parameters: the lattice 
strength, e or Iq, and the norm (number of atoms/total power of the light 
signal), N = J J \^{x,y)\ 2 dxdy. 

Equations JU and J5J represent a class of models that extend to other phys- 
ical media, such as nonlinear photonic crystals and photonic-crystal fibers |17| . 
The results reported below are valid for all these models, i.e., they suggest, 
in particular, the existence of stable vortex complexes of the new types in the 
photonic crystals as well (a vortex with s = 1 was very recently predicted in 
photonic-crystal fibers [T%])- 

We built patterns with vorticity S as sets of M stable "single-cell" solitons 
with s = 0, or compact vortices with s — 1, by imprinting onto the set a phase 
distribution with the phase shift Ac/) — 2wS/M between adjacent solitons. The 
stability of the resulting states was tested in direct simulations. Initial shapes 
of the building blocks with s — and s = 1 were taken as numerical solutions 
of the axially symmetric CPE with an ad hoc potential, U(r) = — ecos(2r) 
for < r < 7r/2, and U(r) = e for r > tt/2, where r = y/x 2 + y 2 , and e is 
the same as in Eq. Q And small random perturbations were further added. 
In cases of stabilities, it was observed that the shapes and structures of the 
vortex and supervortex composed of the building blocks changed hardly in time. 
It implies the vortex and supervortex structures are rather stable, since our 
initial conditions are slightly deviated from the exact stationary solutions. On 
the contrary, the vorticity-carrying phase pattern in unstable states is quickly 
replaced by an evidently chaotic distribution of phase, and shows no trend to 
relax into any simple configuration. The phase instability produces little effect 
on the amplitude shape of the pattern. That is, the profile of \4>\ has hardly 
changed in time. 

Following this procedure, a stable vortex in Eq. (JJ) with S = 1,2,3 and 4 
were created using eight solitons surrounding an empty central cell. We have 
found that a vortex with S — 3 is definitely stable, the one with S = 2 is 
marginally stable. The "marginal stability" features the very slowly (linearly) 
growing perturbation in the phase configuration. Further, vortex solitons were 
also constructed based on 12 single-cell solitons surrounding a cluster of four 
empty cells. In the latter case, the vortices are stable for S — 4, 5, and 6. 
A stable vortex with S = 4 is shown in Figs. H^a) and (b). Monitoring the 
rotation of the local phase vector, which is defined as one with the components 
(Re{ip} lm{ip} / , along a closed path in the counter-clockwise direction 
in Fig. mb), it is easy to see that the corresponding net change of the field's 
phase is indeed Ac/) = 8ir, corresponding to S = 4. 
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Figure 1: (a) A stable vortex soliton with 5 = 4 and N = 33.4; the field \ip(x, y)\ 
is plotted, (b) The field of the local phase vector, (Re{ip} /\iJj\, lm{^p} / \ip\) , at 
centers of the constituent solitons for the vortex complex in (a). In all the cases 
shown here, e = 3. 

The existence and stability of these solitons depend on 5, norm N, and s. 
Keeping e — 3, and varying both N and 5, in the right part of Fig. [3 we 
summarize stability results for the vortex solitons of Eq. Q based on the ring- 
shaped array of 12 single-cell solitons. For the judgement of the stability, we 
have calculated the time evolution of the maximum amplitude of to check 
the collapse, and the time evolution of the root mean square of A<fi — 2irS/12 
for the 12 single-cell solitons to investigate the phase instability, where Acj) is 
the phase difference between adjacent solitons. The numerical simulation was 
performed until t = 150. We have constructed the stability diagram for the 
vortex solitons in Fig. 2 from these two kinds of time evolutions. For all 5, 
the pattern collapses (each individual soliton suffers, as usual, catastrophic self- 
compression) if N exceeds a critical value, N CI = 68.2, that virtually does not 
depend on 5. This can be easily understood, taking into regard that, in the 
present notation, the critical norm for the onset of the collapse of an individual 
soliton in the 2D free-space nonlinear Schrodinger equation is N^ LS ^ w 5.84 
|19| . hence for the set of 12 pulses it predicts N CI « 70.1. The small difference 
from the actual value, N CI — 68.2, is due to the OL. For N < N CI , all the 
patterns with 5 < 3 are unstable, while the one with 5 = 3 is marginal. 

The above results for the vortices based on the ring-shaped arrays of 4, 8 and 
12 single-cell solitons suggest that the stability is not determined by the vorticity 
5 itself, but rather by the phase difference A(f> between adjacent cells. All the 
vortex complexes built as described above are stable for Acf> > ir/2 (including the 
quadrupoles, for which A<p = tt), being marginal for A<p = tt/2. These inferences 
can be easily explained. Indeed, it is well known that sufficiently separated 2D 
solitons with the phase difference A<j> interact repulsively if A(j> > 7r/2, and 
attractively if Atfi < tt/2. On the other hand, it is also known that bound states 
of lattice solitons may be stable only in the case of repulsion [201 - Thus, the 
ring arrays trapped on the lattice have a chance to be stable just in the case of 
A(j) > tt/2. 

It may be relevant to look at several fundamental vortices with 5=1 studied 
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before from the same perspective. Baizakov et al. found the fundamental vortex 
with 5=1, which is arranged as a complex of four major and four minor peaks 
surrounding a central cell of the OL (see Fig. 3(b) in [3]) and they showed that 
the vortex is stable by direct numerical analysis |3J. Alexander et al. found 
numerically a (stable) vortex with 5 = 1 constructed of 4 single cell solitons 
for the PRC model JT], where the phase difference is = ir/2. Yang studied 
the linear stability of the vortex with 5 = 1 and found that the vortex is stable 
in an intermediate range of peak amplitude and it is unstable outside of the 
parameter range |S] . We have found that the vortex with 5=1 constructed of 
4 single cell solitons is marginally stable (very weakly unstable) judging from 
the time evolution of the phase configuration. However, the amplitude profile 
|^| keeps the initial profile stably and the vortex structure is not broken until 
t = 150 for the parameter value of s = 3. In any way, the marginal case of 
A<j) = ir/2 seems to be delicate and might need further investigation. However, 
we do not study it more in detail in this paper, because the main topics are 
higher-order vortex solitons and supervortices. 
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Figure 2: The right side of the diagram shows stability limits for the vortex 
solitons with higher vortices of the type shown in Fig. ^c), for e = 3. Above 
the critical norm (crosses), each soliton in the patterns collapses. The rings are 
stable in the shaded areas below the borders shown by rhombuses. In particular, 
the vortex with 5 = 4 is definitely stable, while the one with 5 = 3 is marginally 
unstable. The left side shows the same limits for densely-packed supervortices 
of the type displayed in Fig. |3 with e = 10. The supervortices are stable in the 
shaded area, which actually includes only 5 = ±1 and 5 = 0. 

Next, we investigate supervortices. To this end, we built a closed array of 12 
compact vortices with s = 1 surrounding an empty cluster of four cells, imprint- 
ing the global vorticity onto it the same way as above (unlike the free space, 
where the vortex solitons is always unstable against azimuthal perturbations, 
in the lattice it may be completely stable, if taken below the collapse threshold 
HUE]). The resulting stationary patterns are, to the best of our knowledge, the 
first example of "supervortex" structures on lattices. On the contrary to the 
vortex rings considered above, the supervortices corresponding to 5 and —5 are 
not equivalent. In the same case, the supervortices with 5 = — 1 and |5| =2 
are stable too, the ones with |5| > 3 are unstable, and the supervortices with 
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Figure 3: A densely-packed stable supervortex with e = 10, S = — 1 and N = 
74.01. (a) The plot of \ip(x,y)\. (b) The area with Re{ip(x,y)} > 0.2 is shaded, 
to highlight the vortex structure of the pattern, corresponding to S = —1. 
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Figure 4: Stable vortex patterns in the photorefractive model (j2J), with Iq = 0.5 
and K = tt/10. (a) A vortex soliton with N = 438.7 and S = 4. (b) A 
supervortex with N = 938.7 and S = —1. The area with Re {ip(x, y)} > 0.1 is 
shaded, to highlight the vortex structures. 



\S\ = 3 are marginal. We have also constructed the simplest supervortices, 
composed of four small vortices which cover four adjacent cells of the square 
lattice (not shown here). It was found that such arrays are stable with S = 0, 
unstable with \S\ =2 (actually, these are vortex quadrupoles) , and marginal 
with S = ±1. 

The stability of the supervortices composed of 12 or 4 individual vortex 
solitons suggest that all the supervortices with the phase difference A0 between 
adjacent individual vortices are unstable in the case of A0 > tt/2, and they 
may be stable if A<fi < tt/2, exactly opposite to the vortex rings considered 
above. This property can be understood too: the sign of the interaction between 
compact vortex solitons differs from that for their zero-spin counterparts by the 
factor of (— l) s |21j . Therefore, contrary to the solitons with s = in the ring 
patterns considered above, the adjacent vortices with s = 1 in the supervortex 
patterns repel each other, provided that A</> < 7r/2, which gives them a chance 
to form stable bound states on the lattice, as per Ref. |20|. 

Supervortices may also exist in a densely-packed form, without the inner 
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hole. As shown in Fig. EJa), it is possible to build such structures, starting with 
a dense set of 16 compact vortex solitons. The phase pattern corresponding 
to the global vorticity S was imprinted, in this case, by setting A0 = irS/2 
in the inner layer (formed of four cells), and Acj) = wS/6 in the outer one, 
that includes 12 cells. Results of the stability investigation for the densely 
packed supervortices are included in Fig. |21 In particular, the collapse threshold 
(above which each individual vortex collapses independently) is nearly constant, 
iV c ( r supcr) « 125. This value can be explained similar to how it was done above 
for the patterns composed of the s = solitons, taking into regard the known 
critical norm for the collapse of an individual vortex soliton with s = 1 |22|. 
The densely packed supervortices with |5 > 2\ are unstable. 

The PRC model gives rise to vortex solitons and supervortices very simi- 
lar to those found in the GPE (JJ, obeying similar phase-stability conditions (of 
course, there is no collapse in the PRC model with the saturable nonlinearity) . 
Figures 4(a) and (b) display a higher-order vortex soliton with S = 4, and a 
supervortex with S = — 1. They are based on 12 single-cell s = or s = 1 
solitons which surround four empty central cells. 

To summarize, we have investigated new types of vortex complexes in a class 
of models that includes attractive BEC in square optical lattices (OLs) and 
photorefractive crystals (PRCs) with the photonic lattice. The same patterns 
are expected in nonlinear photonic crystals and photonic-crystal fibers. The 
complexes include densely packed and ring-shaped higher-order vortices and 
supervortices. They are stable provided that the OL is strong enough. Stability 
limits for the total norm and vorticity of the patterns have been found and 
explained. The experimental creation of these vortex complexes in PRCs is 
quite feasible. 

A larger variety of the complex vortex patterns is expected on triangular and 
hexagonal lattices, and, especially, in the 3D BEC model with the corresponding 
OL, cf. the findings made in the discrete 3D model 22 • A vast "zoo" of 3D 
vortex patterns was very recently reported in Ref. |23|. but they were formed 
of dark-soliton vortex filaments in the repulsive GPE. 

We appreciate discussions with B. B. Baizakov, M. Salerno, and M. Segev. 
The work of B.A.M. was supported, in a part, by the Israel Science Foundation 
through the grant No. 8006/03. 
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